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Q ■ The dynamics of a freely jointed chain in the continuous limit is described by a field 

theory which closely resembles the nonlinear sigma model. The generating functional 

\ ^[J] °f this held theory contains nonholonomic constraints, which are imposed by 

43 ■ 

inserting in the path integral expressing ^f[J] a suitable product of delta functions. 

^ ! 

The same procedure is commonly applied in statistical mechanics in order to enforce 
topological conditions on a system of linked polymers. The disadvantage of this 

^— > . 

method is that the contact with the stochastic process governing the diffusion of the 

C^- 
chain is apparently lost. The main goal of this work is to reestablish this contact. To 

. this purpose, it is shown here that the generating functional ^[J] coincides with the 

£ : 

generating functional of the correlation functions of the solutions of a constrained 
Langevin equation. In the discrete case, this Langevin equation describes as expected 

7— I ■ 

■ the Brownian motion of beads connected together by links of fixed length. 

o : 
o . 
\o . 

m ; INTRODUCTION 

(N ■ 
t— i ; 

The subject of this work is a chain obtained by performing the continuous limit of a 
j> | system of N — 1 links of fixed length a and N beads of constant mass m. In this limit 



^ ■ the number N of beads approaches infinity, the length of the links and the mass of the 

h : 

beads go to zero, while the total length L of the chain remains finite. The dynamics of a 



ater on, mainly the statistical 
[4, 0]. Dynamical models are 



chain with rigid constraints of this kind has been studied in a remarkable series of papers 
ll, l2l, |3| using an approach based on the Langevin equation, 
mechanics of such chains has been investigated, see e. g. 
however interesting by themselves and have also some applications, for instance in modeling 
the response of a chain to mechanical stresses in micromanipulations Q]. 

In Ref. js] the dynamics of the constrained chain has been considered using path integral 



methods. The resulting model is a generalization of the nonlinear sigma model which will 
be called here the generalized nonlinear sigma model or simply GNLcrM. The most striking 



difference from the standard nonlinear sigma model is that in the GIN 



If], 



LcrM the constraint is 



111 ] has been discussed 



]. It has also been shown that it gives the correct equilibrium limit in agreement 
1|. Applications of the GNLcrM have been developed in Refs. Q, Q], computing 



nonholonomic. The relation of the GNLcrM with the Rouse model 
in Ref. 
with Ref. 

for instance the dynamic form factor of the chain in the semiclassical approximation and 
the probability distribution Z(r 12 ) which measures the probability that in a given interval 
of time the average distance between two points of the chain is r±2- 

One point that still needs to be clarified is if the GNLcrM can be related to some stochastic 
process. In fact, the GNLcrM has not been derived starting from a Langevin equation 
and applying for instance the Martin-Siggia-Rose formalism [12| in order to pass the path 
integral formulation. The problem is that this approach becomes cumbersome if one has to 
deal with constraints. For this reason, in [8] the constraints have been added to the path 
integral describing the fluctuations of the beads with the help of an insertion of Dirac delta 
functions. This is a widely exploited procedure in the statistical mechanics of polymers in 
order to impose topological conditions 

To establish a relation between the GNLcrM and a stochastic process is the main goal of 
the present work. To this purpose, after a brief introduction to the GNLcrM, we define in the 
next Section a two dimensional vector field ip v which satisfies a free Langevin equation and 
additional nonholonomic constraints. These are exactly the same constraints which appear 
also in the GLNcrM. Our treatment is limited to two dimensions for simplicity. The gener- 
ating functional ^f[J] of the correlation functions of the fields ip v can be constructed using 
the prescription of [19]. The discretized version of ^f[J] describes the Brownian motion of 
a set of N— beads with diffusion constant D which are connected together by links of fixed 
length. The difference between iff [J] and the generating functional *&[J] of the correlation 
functions of the GNLcrM consists in a functional determinant. We show that this deter- 
minant is trivial by eliminating the constraints using a special set of variables, called here 
pseudo-polar coordinates. As a result we prove the equivalence of ^[J] and &[J] and thus 
the connection of the GNLcrM with a stochastic process of diffusing particles. 
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THE GENERALIZED NONLINEAR SIGMA MODEL AND ITS RELATION TO 

THE LANGEVIN EQUATION 

The starting point of this Section is the generating functional of the GNLSM presented 
in Ref. [8|: 

ijf[j] = J VK(t,s)e- c ^ fdt Io Nds ^ 2 ( t > s )S(\R , (t,s)\ 2 - i) e -/>J>Jf^)^») (i) 

with R = dH/dt and R' = dH/ds. The boundary conditions at t = and t = tf of the 
field R(£, s) are respectively given by R(0, s) = Ro(s) and R(t/,s) = R/(s), where Ro(s) 
and R/(s) represent given static conformations of the chain. For a ring-shaped chain the 
boundary conditions with respect to s are periodic: R(t, s) = R(t, s + L). An open chain 
with fixed ends may be described using instead the boundary conditions: R(£, 0) = rj, 
' \(t, L) = r2, T\ and r2 being the locations of the fixed ends. It was shown in Refs. [8j and 



1 31 ] that the above generating functional describes the dynamics of a closed chain that is 
the continuous version of a freely jointed chain consisting of links and beads. The constant 
c appearing in Eq. (TjQ) is given by: 

M 

Here ks denotes the Boltzmann constant, T is the temperature and r is the relaxation 
time which characterizes the ratio of the decay of the drift velocity of the beads. M and L 
represent the total mass and the total length of the chain respectively. The starting point to 
derive *&[J] is the path integral \I/jv describing the brownian motion of a set of iV particles. 
The rigid constraints, which take into account the fact that these particles form a chain and 
thus are connected together by iV — 1 massless segments of fixed length a, are introduced 
in the path integral with the help of a suitable product of Dirac delta functions. The limit 
of \l/jv from the discrete system to the continuous chain in which iV — > +oo, a — > and 
Na = L can be performed rigorously. The result is the generating functional of Eq. (Tj[|). This 
procedure is different from the usual approach to the dynamics of a chain, which is based on 
a Langevin equation. In this Section we are going to show that the GNLcrM can be related 
to a Langevin equation too. For simplicity, we restrict ourselves to the two dimensional case. 

Since the GNLcrM ignores all interactions, it is natural to suppose that it should be 
related to a Langevin equation with no external forces: 

<p„ = v (3) 
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where <*p v is a two dimensional vector field and v is a white noise source, whose components 
i/W, i — 1,2 satisfy the basic correlation functions: 

<!/»(*,«))> = (4) 

(i/«(f, (f.aO) = ^(t-O^a-a') i,j = 1,2 (5) 

One may also expect that, together with Eq. (j3J), the field ^ must also satisfy the constraint: 

tf, 2 = 1 (6) 
The generating functional *&[J] of the correlation functions of the field y> v is then given by 

The meaning of the statistical sum in the right hand side of the above equation becomes 
clear if we rewrite it as follows: 

f[J]= I Vv I VRe- c ti f * ft*"* SiR-tpJefi* (8) 

J JR' 2 =1 

The path integration over u is now unconstrained, while that over the new field R is limited 
to the configurations which are of the form: 

R(t, s) = / du(cos(j)(tj u), sin <p(t, u)) + Ro(t) (9) 

where Ro(i) is independent of s. The only left degree of freedom is the angle 4>(t, s). 

The generating functional &[J] of Eq. ([T]) differs from ^f[J] due to the presence of the 
functional Dirac delta function 5(R' 2 — 1). As a matter of fact, it is easy to show that: 

V[J] = J Vv J DRe- c ^*^^5(R' 2 - 1)5(R- v )e~ & * $ (10) 

The connection with the Langevin equation ([3]) is made by noticing that, for any solution 
«*p v of that equation it is possible to write the formula: 

8(R - u) = detT^^R - if v ) (11) 

Applying Eq. (TTTT) to Eq. ffTUj) we obtain, up to an irrelevant constant: 

*[J] = /" DRPi/e'^o^'^^^R - <p„)5(R' 2 - i) e -fo f dtf L ds3.n ^ 
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As already announced, this expression of the generating functional ^[J] differs from iff [J] 
only by the fact that the condition R' 2 = 1 is imposed with the help of the delta function 
<5(R' 2 — 1). In the next Sections the degrees of freedom which are frozen by the condition 
R' 2 = 1 will be projected out from the path integral ( 1T21) and it will be shown that what re- 
mains is exactly the generating functional ^f[J] related to the constrained stochastic process 
of Eqs. © and ©. 



THE DISCRETE GENERATING FUNCTIONAL IN PSEUDO POLAR 

COORDINATES 

As a first step to show the equivalence of the generating functionals *ff[J] and ^f[J] we 
replace the continuous variables s and t with discrete variables s m and t n , with < m < M 
and < n < N. The spacings in the discrete s and t— lines are respectively given by: 

s m ~ s-m-i = a m = 2,...,M (13) 

tn-tn-i = b n = 2,...,N (14) 

where a and b are supposed to be very small. The continuous limit is recovered in the limit 
M, N — > +oo, a, b — > and Ma = L, Nb = tf. To simplify formulas, it will be used in 
the following the shorthand notation: 

In this way the discrete version of the constraint R' 2 (t,s) = 1 is replaced by the set of 
conditions: 



[Rum — Rn(m-l)) _ 71 — 1, . . . , iV 

^2 



(16) 



a m = 2, . . . , M 

With the above settings the generating functional \P[J] of Eq. f fT2l) may be rewritten as 
follows 0: 

-+oo 



v&[Jl = lim lim / 

N^ooM^oo / 



r, 



exp 



-abc v, 



2 

nrn 



X 



l[5(R 

nm ^Pi>,nm) 6Xp \ db ^ ^ J ram R nm / 
n,m (. n,m J 



2 _ / Rnm R' 



X 

a \ a 

n m=2 
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Let us also note in the last line of equation ( TT7I) the normalization factor Y[ n Ylm=2 a ^ n ^ ne 
definition of the delta function imposing the constraints. The reason of this normalization 
will be clear later. Without the constraints, the above equation would describe a discrete 
chain of N — 1 segments of length a and N beads of mass m respectively which perform a 
Brownian motion. The diffusion constant D is recovered from the parameter c of Eq. (J2]) 
as follows. First of all, we note that ca = ik ^ Tr ^a. The ratio ^ is nothing but the linear 
density of mass along the chain, so that is equal to the mass m of a single bead, i. e.: 
= m. As a consequence, ca = 4k ™ Tr . At this point we remember that the mobility of a 
particle /i may be expressed in terms of the particle mass m and of the relaxation time r as 
follows: ™ = -. Thus, ca = 4k ^ Tfl ■ Due to the fact that D = k B Tfi, it is possible to write 

ca = Id- 

To eliminate the constraints (fl6l) . we pass to a new set of coordinates, which in the 
following will be called pseudo-polar: 

M 

R-nm = ^2 lnm ' ( COS ^m',Sm0 nm /) (18) 
m'=l 

The ranges of variation of the variables l mn and (p nm are respectively given by: 

< l mn < +00 < (j) nm < 27T (19) 

The coordinate l nm for n = 1, . . . ,N and m = 2, . . . , M, describes the length of the m— th 
segment at the instant t n . The coordinate l n i is very special, because it gives the position 
with respect to the origin of the reference system of the first bead in the chain at the time 
t n . Finally, the angles <p nm tell us how the N — 1 segments are reciprocally oriented. After 
the transformation (fTBl . the vector R nm depends on the variables l nm and nm , i. e.: 

Rnm Rnm({Lm}] ^nli {0nm}) (^0) 

where {l nm } is the set of all / nm 's for which m ^ 2 and {0 nm } is the set of all nm 's. 
Analogously, we denote with {R nm } the set of all R n m's for m = 1, . . . , M and n = 1, . . . , N. 
We are now able to explain the reason of the normalization factor Y[ n Ylm=2 a ^ n 03) • m 
the pseudo-polar variables the constraints ([161) become: = 1. The factor - is necessary 
in order to normalize the delta functions imposing these constraints. As a matter of fact, it 
is possible to check that | / + °° dl nm 5{ lj: ^- — 1) = 1. 
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In order to perform the transformations (ITS]) in the expression of the generating functional 
^[J] of Eq. f|T7|) . we need to compute the associated Jacobian determinant. In the rest of 
this Section we will prove for a general functional /({R nm }) the following formula: 



[J c/R nm /({R nm }) 



+00 



2tt 



} J <^nm / Y\_ dfinmf ({Rnm({lnm} , l-ril] {finm})}) JnM 



where the Jacobian Jnm of the transformation (|T8|) is given by: 

JNM{\)nm\i ^nli {4>nm}) = \\ Uihdn(M-l) ' ' ' Inl 



(21) 



(22) 



Let's show that Jnm is really that given in Eq. ff22l . In order to proceed, it is convenient 
to introduce the components Xnm and Xnm of the vectors R nm , i. e. R nm = (xnm, Xnm) ■ 
Thus, Eq. f fT8|) becomes: 



in 



x 



(i) 



/ COS( 



X 



(2) 



m'=l 

and Jtvm m ay be written as follows: 

JNM({lnm}, Inl', {4>nm}) = det 



m'=l 



(23) 



dx (1) 


dx {2) 






dx W 


dx {2) 


d <t> n 'm> 


9( t> n 'm' 



Strictly speaking, Jnm is the determinant of a block matrix A r , 



(24) 



with composite 



indices nm and n'm'. A nm . n i m i 



is composed by four NM x NM matrices, since n, n' 



a (*) 

1.....N and m, m! = 1.....M. Due to the fact that Ip^ 



dx (i) 



for i = 1,2 if 



n ytz n', A nm . n r m i is a block diagonal matrix. As a consequence, it is possible to write its 
determinant as follows: 

Jnm = Y\_ J' nM (25) 



where 



J, 



nM 



det 



dx W 


dx (2) 


9l nm> 


dl n m > 


dx W 


dx {2) 


d( t>nm> 


d( t> n m> 



(26) 



Using Eqs. (1231 . one finds after a few calculations that J„m is the determinant of the block 
matrix: 

A(n) B{n) 
C(n) D(n) 



J, 



nM 



det 



(27) 
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A(n), B(n), C(n), D(n) are lower triangular M x M matrices with elements: 



A m m'{j^) @mm' COS (finm' 
Cmm' (^) Inm' 9 mm' sin (j) nm ' 



Bmm' \n) = 9 mm > sin <\> nm i 

D m m' (jl) Inm' 9 mm' COS (f)nm' 



(28) 
(29) 



Here the matrix 6 mm > denotes the discrete equivalent of the Heaviside theta-function: 



Omm' = 1 if m! < m 



Omm' = if m' > m 



(30) 
(31) 



If the matrices A(n), B(n),C(n), D(n) would commute, one could use a known theorem of 
linear algebra and write: J n u = det(A(n)D(n) — B(n)C(n)). In our case these matrices do 
not commute, but it is still possible to compute the determinant J n u by induction on M. 
If M — 1 it is easy to show that: 

■Jnl — ^nl (32) 

Next, we prove that 

JnM = InMJn(M-l) (33) 

To this purpose, it will be convenient to introduce new indices a, (3 = 1, . . . , M — 1. At this 
point, we note that the M— th column of the 2M x 2M block matrix whose determinant we 
wish to compute in Eq. ( 1271) has only two elements which are not zero. Thus, we expand J n u 
with respect to the M— th column. Taking into account the necessary permutations and the 
fact that the two nonvanishing elements are A MM {n) = cos 0„m and Cmm(«) = —l nm sin 4>nM 
we obtain: 

aj3 cos n/3 d a/3 sin (j) nf3 
JnM = cos0„ A/ det -l n p6 a/3 sin <f) nf3 l n p6 a p cos <j) n p 

— InfiQMp Sin 4> n Infi^Mf! COS (f) n p InM COS (j) n M 

9 a p cos <j) n p 6 a/3 sin (j) nf3 
< i; V /det M p cos (j) n j3 9m p sin n/3 sm<fi nM (34) 
-lnp9 a psva.(j) n p l n p9 aP cos 4> n p 

The determinants of the remaining two (2M — 1) x (2M — 1) matrices may be expanded 
according to the (2M — 1)— th column, because these columns contain only one nonvanishing 



+ 



1) M Lm sin 
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element. After simple calculations one finds: 



J, 



nM — l nM 



l„ M det 



which is exactly Eq. (1331) because 



J, 



n(M-l) 



det 



a p cos (p nf 3 9 a/ 3 sin <\> nfi 

-InpOap Sin (j) n p InpQap COS n/3 



a/3 cos 4> n p 9 a(S sin n/3 

-InfiQuP Sin n( g InfiOap COS n/ 3 



(35) 



Using Eqs. ( 1321) and ( 1331) it is easy to show by induction that J u m = l n M^n{M-i) • 
a straightforward application of Eq. (1251) it is now possible to prove Eq. ( 1221) . 



(36) 



Z n i- With 



RECOVERING THE GENERATING FUNCTIONAL *[J] OF THE 
CONSTRAINED STOCHASTIC PROCESS OF EQS. (Sj) ([6]) 

Let's now go back to the generating functional ^[J] of Eq. (|T7|) . After the change 
of variables (fT8l) . the delta functions imposing the constraints simplify as follows: 
^ ^ |Rnm R„( m -i)l _ ]_ j _ § ( hoi — \\ _ Further simplifications are obtained after apply- 



/;2 

ing the two delta function identities o ( 



- 1 



a 2 S(l 



2 

nm 



and 



2~ [S(l n m — a ) + 5(Z nm + a )]- Remembering that in our case l nm > 0, it is possible to put: 
<5(Z^ m — a 2 ) = ^5(l nm — a). As a consequence, the expression of the generating functional 
*&[J] in pseudo-polar coordinates becomes: 



^[J] = lim lim 

N—>oo M^oc 



2tt 



JJ dv nm / Yl dl nm / JJ #nm exp < -abc ^ v l 

J n,m ^ n,m ^ n,m K. n,m 

J~J 5(R. nm ("[/j lm }, l n li {0nm}) ^Pis,nm) GXp \ flfe ^ ^ Jnm * Rnm({^nm}) ^nli {0nm}) 



X I I i nl 

II 



M 



nru« 



— a) 



n m=2 



M 



m=2 



nM 



ill 



(37) 



In writing the above equation we have separated from the Jacobian determinant Jnm the 
contribution coming from the l' n iS, because these quantities denote the positions with respect 
to the origin of the first bead at different times t n 's and are thus not fixed by the constraints. 
The integration in Eq. (1371) over the Z nm 's, for n = 1, . . . , N and m — 2, . . . , M, produces as 
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a result: 



IT dv nm / TT Mrtiinl / IT d(j) nm exp I -dbc V I/* r 
n,m u n u n,m ^ n,m 

x ]^[5(R nm ({a},/ n i; {<p nm }) 

n,m ^ n,m J 

Here the symbol {a} denotes the set of all Z nm 's for m ^ 2 after the imposition of the 
constraints Z nm = a. We can now rewrite Eq. (1381) as an integral over a restricted domain D: 

*[J] = lim lim / TTd/ nm d^ nm / TT cfe/ nm exp <| -a&c V" i% > TT / nM ••• Li 

n,m n,m ^ n,m J n 

X JT 5(R, nm ({/ rim }, {4>nm}) ^Pu,nm) GXp < 0,6 ^ ^ Jnm ' Rnm ( { ^nm } > ^nli {^ram}) ( (39) 



where D is the domain of all i nm 's and nm 's with the constraints Z nm = a for m — 2, 



M 



and n = 1, . . . , N: 



D = ( 



1 1 



m — 2, .... M and n — 1, 



AT 



M and n = 1, . . . , JV 



< i n i < +oo n = l,...,iV 

< <Pnrn < 2tT m = 1, . 

At this point, using Eqs. (I2lj) and (|22l) we go back to cartesian coordinates: 

r r r+oo ^ ^ f 

*[J] = lim lim / T\dR nm l TT du nm exp <^ -a&c i/^„ 

N-^oo M— >oo If) J —no 

n,m n.m \ n,rn 

x [jd(R nm - ¥V, nm ) exp < a& ^ Jnm ■ Rnm f 

n,m ^ n,m J 



(40) 



(41) 



The domain D in cartesian coordinates is given by all R nm 's in the two dimensional plane 
subjected to the constraints (fTBj) : 



Wi.m~R-n(m-l) I 



m= 1,...,M n= l,...,iV 

1 m = 2 M n = l A" 



(42) 



Finally, we rewrite the path integral in Eq. (j4"Ii) in its continuous form. The result is: 

Jr' 2 =i </ 



(43) 



The right hand side of the above equation coincides exactly with the right hand side of 
Eq. (jSJ). This proves the equivalence between the generating functional \P[J] of the GNLcrM 
and the generating functional St [J] of the stochastic process of Eqs. ©-([6]). 
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(44) 



CONCLUSIONS 

In this work it has been shown that the GNLcrM is related to a stochastic process which, 
after discretization, describes the Brownian motion of N beads subjected to the constraints 
(fj~6l) . These constraints enforce the conditions that the links connecting the beads are of 
fixed length. More in details, it has been proved that the generating functional *&[J] of 
the GNLcrM coincides with the generating functional iff[J] of the solutions of the Langevin 
equation ([3]) and of the constraint ([6]). The fact that the two functionals are equal was not a 
priori obvious, because they differ by the delta function S(R' 2 — 1) which contains quadratic 
powers of the fields. If 8(g(R)) is a delta function imposing the condition g(R) = 0, then in 
general the following identity is valid: 

5g{R) 

lssxj i^rtjuei 

lg(R)=0 

If in our case the functional determinant appearing in the right hand side of Eq. (jUJ) would 
be not trivial, then there would be no chance that (|S]) and (1121) coincide. Luckily, it turns 
out that, after passing to the pseudo-polar coordinates (TTg|) . the delta function <5(R' 2 — 1) 
produces just a functional determinant which is a trivial constant. 
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